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ABSTRACT 

We present a new technique for obtaining model fittings to very long baseline interferometric images 
of astrophysical jets. The method minimizes a performance function proportional to the sum of the 
squared difference between the model and observed images. The model image is constructed by summing 
N s elliptical Gaussian sources characterized by six parameters: two-dimensional peak position, peak 
intensity, eccentricity, amplitude and orientation angle of the major axis. We present results for the fitting 
of two main benchmark jets: the first, constructed from three individual Gaussian sources, the second 
formed by five Gaussian sources. Both jets were analyzed by our cross-entropy technique in finite and 
infinite signal-to-noise regimes, the background noise chosen to mimic that found in interferometric radio 
maps. Those images were constructed to simulate most of the conditions encountered in interferometric 
images of active galactic nuclei. We show that the cross-entropy technique is capable of recovering the 
parameters of the sources with a similar accuracy to that obtained from the very traditional Astronomical 
Image Processing System Package task IMFIT when the image is relatively simple (e.g., few components). 
For more complex interferometric maps, our method displays superior performance in recovering the 
parameters of the jet components. Our methodology is also able to show quantitatively the number of 
individual components present in an image. An additional application of the cross-entropy technique 
to a real image of a BL Lac object is shown and discussed. Our results indicate that our cross-entropy 
model-fitting technique must be used in situations involving the analysis of complex emission regions 
having more than three sources, even though it is substantially slower than current model fitting tasks 
(at least 10,000 times slower for a single processor, depending on the number of sources to be optimized). 
As in the case of any model fitting performed in the image plane, caution is required in analyzing images 
constructed from a poorly sampled (u, v) plane. 

Subject headings: galaxies: jets — ISM: jets and outflows — methods: data analysis — methods: 
numerical — methods: statistical — techniques: interferometric 



1. INTRODUCTION 

Very long baseline interferometry (VLBI) is to date 
the only tool available for the direct study of the struc- 
ture and evolution of active galactic nuclei (AGNs) on 
parsec scales. While early observations were restricted 
to a few radio telescopes and required model fitting 
of th e visibility functio n to obtain the source struc- 
ture dCohen et all fl97Th . good coverage of the (u, v) 
plane of modern arrays and the development of efficient 
imaging techniques guarantee excellent maps of these 



sources ([Rogers et~al"1ll974 iReadhead fe Wilkinson! 119781: 
ICornwellloH TP carson fc Readheaof L984: She pherd et all 
Il994ft . 

Data with sparse (u, v) plane coverage or calibra- 
tion problems can produce images inappropriate for the 
study of quantitative aspects of AGNs, such as posi- 
tions and proper motions of jet features (e.g., iPearsonl 
1999;). In this case, the analysis should be made by 
fitting a discrete number of model components (gener- 
ally two-dimensiona l Gaussian) to the vis ibility data (e.g., 
ICarrara et all Il993t iKovalev et all 120051 ) , frequently cor- 
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recte d with self-calibrat ed techniques in the imaging pro- 
cess (je.g.. Person! 1 19951 ). Fits in the (u,v) plane have 
some advantages, such as error recognition, lack of spu- 
rious CLEA N artifacts, and full angular reso lution capa- 
bility (e.g., WeMS^MM^^^MlMM)- However, 
(u, v) components do not necessarily correspond to real 
features in the image, but may simply be mathematical 
artifacts need ed to reproduce complex brightness visibi l- 
ity data (e.g.. lKellermann et al.ll2004t iLister et al.ll2009bD . 
Besides, errors detected in the complex visibility data will 
not always p roduce measurable effects on the image plane 
(|EkersHl999f ). which means that analyses performed on the 
image plane are not necessarily less reliable than those car- 
ried out on the (u, v) plane. Despite the intrinsic advan- 
tages and disadvantages, many works have used the image 
plane to extract information concerning the structural pa- 
rameters of the observed sources (e.g., Kellermann et al.l 
l2004HWrobel et al.ll2004t iRevnolds et ai.ll2009ft . 

Regardless of the plane where the analysis is taking 
place, the model fitting parameters are usually estimated 
from some maximum likelihood estimator, which finds the 
best set of model parameters that minimize the residual 
differences between synthetic and real data. The conver- 
gence of such algorithms usually depends strongly on the 
initial estimation of the parameters and they are prone 
to finding non-global minimum solutions, especially if the 
object to be modeled is complex. In addition, the model- 
fitting routines generally used in interferometric astron- 
omy (e.g., IMFIT, JMFIT, etc.) suffer from a limitation 
concerning the number of components to be fitted simul- 
taneously in the image. 

To overcome those limitations, we introduce a new and 
powerful technique: the cross-entropy method for global 
continuous optimization (hereafter CE). The method uses 
the image as input data, and searches for the optimal 
model parameters, selecting the best candidates among 
all solutions generated in each iteration and constructing 
new ones from them. There is no limitation to the number 
of Gaussian model-fitting components, and the algorithm 
is able to place self-consistently a practical limit on the 
number of model-fitting components, as will be discussed 
later in this work. 

CE analysis was originally used in the optimization 
of complex comp uter simulation m odels involving rare 
events simu l ations ([Rubinstein!! 1 9971 ) . and was modified by 
lRubinste"inl (|1999D to deal with continuous multi-extremal 
and discrete combinatorial optimization problems. Its the- 
or etical asymptoti c conve rgence has b e en de monstrated 
bv lMargolinl ()2004D . while IKroese et al.l (|2006l) studied its 
efficiency in sol ving continuous mult i-extremal optimiza- 
tion problems. ICaproni et al.l (|2009f ) successfully applied 
the CE technique to determi ne precession m o del pa rame- 
ters of relativistic jets, while iMonteiro et al.l (|201dl ) stud- 
ied Galactic open clusters from CE optimization in color- 
magnitude diagrams. Other examples of applicat ion of the 
CE method are provided in lde Boer et all (|2005l ). 

The basic procedures inv olved CE optimizat ion can be 
summarized as follows (e.g.. IKroese et aLll2006l) : 

1. Random generation of the initial parameter sample, 
obeying pre-defined criteria; 

2. Selection of the best samples based on some math- 



ematical criterion; 

3. Random generation of updated parameter samples 
from the previous best candidates to be evaluated 
in the next iteration; 

4. Optimization process repeats steps (2) and (3) until 
a pre-specified stopping criterion is fulfilled. 

In this work, we validate our CE model-fitting algorithm 
from a variety of benchmark tests built from synthetic jet 
components and apply it to a real image. The paper is 
structured as follows: in Section 2, we introduce the CE 
algorithm and its application to the problem of modeling 
interferometric radio images of astrophysical jets. The val- 
idation tests and their respective optimization results are 
discussed in Section 3. The application of our technique 
to a real image is presented in Section 4. Conclusions are 
presented in Section 5. 

2. THE CROSS-ENTROPY METHOD FOR CONTINUOUS 
OPTIMIZATION 

In this section, we introduce the cross-entropy method, 
describing how to estimate the model-fitting parameters 
of jet knots present in the radio interferometric maps of 
astrophysical sources. 

2.1. Global View of the Optimization Process 

Let us consider an interferometric image composed of 
N x x N y pixels, where N x and N y are, respectively, the 
number of pixels in right ascension and declination coor- 
dinates. Assuming that such data can be described by 
an analytical model with JV p parameters p\ , pi , ...,pat p , we 
can use the CE optimization method to find the set of 
parameters x* = (pl,P2, ■■■iP*n p ) for which the model pro- 
vides the best des cription of the data (Rubinstein 1999; 
IKroese et aHl2006l ). 

The optimization process is performed by randomly 
building N independent sets of model parameters X = 
(xi,x 2 , ...,Xjv), where Xj = (pu,P2i, — ,PNpi), and min- 
imizing a performance function S'(x) used to check the 
quality of the fit during the run process. In an ideal situ- 
ation in which all parameters have converged to the exact 
solution, we must obtain S^x*) — > 0. 

To find the best solution from CE optimization, we start 
by defining the parameter range in which the algorithm 
will search for the best candidates: p™ n < Pj(k) < p™ ax , 
where k represents the iteration number. Introducing 
PM = {Pf +Pj max )/2 and Oj-(0) = (pf™ - pf n )/2, 
we can compute X(0) from: 



X ij (0)=p j {0) + a j {0)G i 



(1) 



where Cry is an N x N p matrix with random numbers 
generated from a zero-mean normal distribution with stan- 
dard deviation of unity. 

The next step is to calculate Si(Q) for each set of Xj(0), 
ordering them according to increasing values of Si. Then 
the first iV e iit e set of parameters is selected, i.e. the iVeiite- 
sample with the lowest S'-values, which will be labeled as 
the elite sample matrix X olltc (0). 

We then determine the mean and standard deviation of 
the elite sample, pf lte (0) and af lte (0) respectively, as: 



VLBI images and cross-entropy optimizer 



3 



pf te (0) 



lite 



iV e ute 

»=i 



nit 



(0), 



(2) 



(Tf tC (0) 



1 Welrt, 



pf te (0)] : 



The matrix X at the next iteration is determined as: 



(3) 



^(l)=pf te (0) + 



elite 



(0)G. 



o > 



(4) 



This process is repeated from Equation (2), with Gij re- 
generated at each iteration. The optimization stops when 
the maximum number of iterations k max is reached. 

In order to prevent convergence to a sub-optimal so- 
lution du e to the intr in sic ra pid convergence of the CE 
method, iKroese et all (|2006ft suggested the implemen- 
tation of a fixed smoothing scheme for p° lltc ' s (fc) and 



elite, s 



(ft): 



pf te ' B (k) 



apf tc (fc) + (l-a)pf te (fc-l), (5) 



af tc ' s {k) = a d (k)af tc (k) + [1 - a d (k)} af tc (k - 1), (6) 

where a is a smoothing constant parameter (0 < a < 1) 
and ad(fc) is a dynamic smoothing parameter at the fcth 
iteration: 



Oia(k) = a — a (l — k 1 ) , 



(7) 



where q is an in teger typically between 5 and 10 
(jKroese et al.ll2006ft . 

2.2. Defining the Performance Function 

In the optimization problem we need to define a perfor- 
mance function based on the desired characteristics of the 
solution. Usually for continuous problems this is done by 
defining a likelihood function and then maximizing it, or 
requiring that the sum of the residuals squared be mini- 
mal. 

In this work, we adopt a combination of the sum of the 
squared residuals and their respective variance. The mo- 
tivation for this procedure is that there is no parametric 
form for determining the error for a given pixel of the inter- 
ferometric radio images as is the case for typical photon- 
counting errors which are known to be Poissonian. There- 
fore we assumed that a good fit of jet sources would essen- 
tially produce small residuals, as well as a uniform spatial 
distribution of them, thus with minimal variance. 

Let the quadratic residual R m (k) at a given pixel m and 
iteration k be defined as the squared difference between the 
observed image I m and the generated model image M m (k) 
at a fc-iteration, i.e. R m (k) = [I m — M m (k)] 2 . The mean 
square residual value of the model fitting R(k) can be cal- 
culated from: 



R{k) 



pixel 



Npixcl 



m—1 



(8) 



We choose to rank our tentative model images obtained 
from the 6N S parameters Xj(fe) at iteration k through the 
performance function: 



S P rod(xi,k) = R(k) x 



pixel 



(R m (k)-R(k)Y 



m—1 



(9) 

in which the product between R and the variance-like fac- 
tor was used. Note also that, if I m and M m are expressed 
in terms of Jy beam -1 , S'prod has units of Jy 6 beam . 

It is important to emphasize that we have tested other 
functional forms for the CE performance function but 
Equation (9) was more efficient in more complex problems. 
We believe that the difference in performance is essentially 
due to the fact that Equation (9) transmits directly any 
change in the mean and variance of the residuals in all 
iteration steps of the optimization process. 

2.3. Estimation of the model fitting parameters and their 

uncertainties 

As commented previously, CE optimization generates 
TV random tentative solutions at each iteration k, select- 
ing the best N c \n c set of model parameters in terms of the 
values of S pro d- 

In all optimization process presented in this work, wc 
first chose the number of sources N s to be model fitted 
on the image and applied CE optimization three times for 
each choice of N s . After finalizing the three optimization 
runs, we selected the optimization that best minimized 
<Sprod in order to determine the best values of the 6iV s 
Gaussian parameters p* , as well as their respective uncer- 
tainties <7 p * (i varying from 1 to 6N S ) as follows: 



Pi = Y WkPlk Wk 



(10) 



fc=i 



, k=l 



and 



fcmax /femax \ 

<r = X>*(Pi*-K) 2 X>* . ( n ) 

k=l \ k=l ) 

where fc max is the maximum number of iterations used 
in the optimized model fitting, and pik represents the set 
of model parameters that produce the minimum value of 
■S'prod,, among all tentative solutions at iteration k and 
Wk = Sp TO< i ■ The power index -2 in the definition of Wk 
was adopted in order to make the tentative solutions with 
the lowest values of S'prod, more important in the calcula- 
tion of p* . 

The same procedure was used to estimate the value of 
S'prod associated with n*, except for Wk — cr^ 4 , where 

<7q is the variance of SV^e-d among all tentative solu- 
tions at an iteration k. 

3. VALIDATING THE CE MODEL-FITTING METHOD 

To validate the fitting technique, we used synthetic im- 
ages generated from a previously known set of parameters. 
Such images were composed of N s two-dimensional ellipti- 
cal Gaussian sources, characterized by the Gaussian center 
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peak positions (x, y), semi- major axis a, eccentricity e, po- 
sition angle 9, and peak intensity I . 

We created two synthetic jets, labeled Jl and J2, with a 
bright main central source (emulating an unresolved core 
as seen in real interferomctric images) and additional com- 
ponents with a range of sizes and with different orienta- 
tions, intensities, and degrees of superposition. We took 
them as being representative of typical jet sources found 
in the literature. 

We also ran simpler tests, such as single source, fully 
separated set of sources, etc., but we do not present them 
here because they do not give any more insight into the 
performance of the method. In all these simple cases the 
method performed extremely well, finding the correct so- 
lution with considerable ease. 

Our method was implemented in such a way that the 
fitting was entirely automatic. This has opened the pos- 
sibility of performing fits for a wide range and number of 
components, analyzing their quality a posteriori based on 
more quantifiable characteristics. In particular, we have 
used mainly the best value of the fitness function S* rod (xi) 
and the characteristics of the residual maps to decide the 
ideal number of sources, as can be seen in the following 
sections. 

It is important to emphasize that we have adopted 
fc max = 200CW S , a = 0.9, q = 5 N = N (6N S ) 2 (with 
30 < No < 50) and N e i ite = 0.05AQ in all validation tests 
presented in this work, in agreement with those suggested 
bv lKroese et~aTl (|2006f) and Caproni et al. (2009). Addi- 
tional tests were made using different values for those CE 
parameters but they did not improve the overall perfor- 
mance of our technique. The CE optimizations presented 
in next sections were run on a 2.7 GHz processor, taking 
typically from a couple of minutes up to days to be fully 
completed, depending on the number of elliptical Gaussian 
components assumed in the model fitting. 

3.1. Noiseless Synthetic Benchmark Tests 

This section is devoted to showing the performance of 
our CE model fitting in ideal situations where there is 
no background noise present, i.e., the signal-to-noise ratio 
(S/N) tends to infinity. 

3.1.1. The Noiseless Synthetic Jets Jl and J2 

Jet Jl contains three elliptical Gaussian sources (-/V s = 
3) for which peak intensity /size decreases/increases with 
the core-component distance. Their respective parameters 
are shown in Table [TJ The synthetic image of Jl with a 
size of 30 x 90 pixels in right ascension and declination, 
respectively, is displayed in Figure [T] While the brightest 
component (mimicking the core of the object) is relatively 
detached from the others, the remaining two components 
present a moderate degree of superposition, mainly be- 
cause of the size of component 3. 

Jet J2 is composed of five distinct elliptical Gaussian 
sources (N s = 5), with a small central source and the 
others at increasing distances, and larger sizes, as usually 
observed in astrophysical jets. The source parameters are 
listed in Table [5J while the resulting map covering an area 
of 50 x 60 pixels in right-ascension-declination plane is 
shown in Figure [2j The distances between the compo- 

1 ICaproni et all <201ll) adopted fc max = 5000 and N = 50(6JV S ) 2 . 



nents were chosen to introduce different degrees of spatial 
superposition. We also put a weaker source between two 
bright ones (source number 3) with high degree of super- 
position to evaluate the ability of the method to separate 
sources with intensities slightly above the wing intensities 
of neighboring sources. 

3.1.2. The CE Model- fitting Results for Jl and J 2 
Noiseless Maps 

The CE model-fitting process generates elliptical Gaus- 
sian sources with peak positions within the observed im- 
age boundaries and from a large span of parameter space 
as follows: 2.5 < a(pixel) < 40, < e < 0.95 and 
—90° < 9 < 90° (positive values for north to east di- 
rection) for both Jl and J2 jets. In relation to the peak 
intensity, the limits are 0.0002 < Iq (Jy beam _1 )< 2.1 for 
Jl and 0.0001 < I Q (Jy beam" 1 )^ 1.2 for J2. 

The values of the model fitting parameters obtained 
from our CE global optimization technique for Jl and J2 
benchmark tests are shown, respectively, in Tables [3] and 
01 Those values, as well as their respective uncertainties, 
were calculated from the procedures described in Section 
2.3. The number of sources adopted in each optimization 
was varied from two to seven for Jl and three to seven 
for J2 (justification for this procedure is given in the next 
subsection). 

A quick comparison between Tables [5] and U] reveals 
that our algorithm fits tentatively the most intense sources 
present in the image if N s < 5. In other words, compo- 
nents 1 (core) , 2 and 4 are always present in the modeling, 
even though their recovered parameters are not necessarily 
correct in all cases. If we employ iV s > 5, the algorithm 
places the extra Gaussian sources close to the location of 
the most intense sources. Except for the threshold value 
N s = 3, the same considerations are valid for the compar- 
ison between Tables Q] and [3] related to jet Jl. 

Comparison between the values of Table [5] and the five- 
source entries of Table @] reveals an excellent agreement 
between them. The same is true in the case of jet Jl 
(compare the three-source entries of Tables Q] and [3]) , for 
which our algorithm recovered the structural parameters 
of the three sources with a tantalizing zero relative error. 

3.1.3. Estimating the Number of Jet Components from 
the CE Model-Fitting 

There is an important aspect that must be addressed 
when any model fitting technique is applied to an interfer- 
omctric image: how many components are supposed to be 
used in the fit? Unfortunately the answer is that there is 
no ultimate rule to apply in order to determine the num- 
ber of components to be fitted in the image (even a mini- 
mum or a maximum limit). Perhaps the most reasonable 
and conservative criterion could be the one based on the 
Occam's Razor principle, adopting the lowest number of 
sources for a reasonable fit. This strategy can be found in 
some previo us works that deals w ith model fitting of VLBI 
images (e.g., iHoman et aT1l2001l) . 

To verify if our algorithm can determine quantitatively 
the number of sources, or at least the minimum number of 
components that should be assumed in the modeling, we 
have varied the number of Gaussian sources in the CE opti- 
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mization procedure. We show in the top panels of Figures 
[3] and [4] the one-sixth valuetQ of the performance function 
S* rod as a function of the number of sources assumed, re- 
spectively, in the model fitting of Jl and 32. The error 
bars are large and only allow us to conclude that there 
is possibly a plateau- like structure after N s ~ 3 for J 1 
and N s ~ 5 for J2. This suggests a crude estimate for 
the minimum value for the number of sources present in 
the image. This is a slight advance in relation to previous 
works since our technique allowed us to put a lower limit 
on the number of the sources, based only on the values of 
our performance function. 

With the aim of improving our estimate of N s , we plot- 
ted in the bottom panels of Figures [3] and |4] the mean and 
the maximum values of the residual images obtained from 
our CE model fittings for Jl and J 2 respectively. These 
quantities indicate that there is a minimum at N s = 3 for 
Jl and N s = 5 and 7 for J2, which means we could not 
unequivocally constrain N s for jet J2. However, looking 
carefully at Table @J we can see that the same Gaussian 
sources for N a = 5 were also found in the model fitting 
employing N s — 7. The extra two components are too 
dim, compatible with zero-intensity taking into account 
the uncertainties. Therefore, the ambiguity in the correct 
number of the components is definitely eliminated. 

It is important to emphasize that similar behavior was 
found for jet Jl, the three sources for N s = 3 having also 
been detected by our CE optimizations in the cases of 
N a = 4, 5 and 6. Thus, we believe that the procedures 
mentioned above can be very useful and effective in putting 
a quantitative constraint on the number of sources present 
in any interferometric image. 

3.2. Noise Validation Test 

In the previous section, we checked the performance of 
our CE model-fitting technique in recovering the struc- 
tural parameters of synthetic sources in the limit of infi- 
nite S/N. However, this idealized situation is definitely far 
from those found when dealing with real interferometric 
radio maps, in which the S/R has a finite value. Thus, it 
is necessary to verify if our technique also works in such 
conditions. 

The question that arises is how to perform tests that 
include a realistic noise background. In the VLBI images 
obtained by the CLEAN method the presence of struc- 
tured noise is evident, in other words, noise not originat- 
ing in Gaussian or Poissonian distributions. Because there 
is no simple parametric expression for the noise in those 
images, we carried out the following procedure to obtain 
our finite-S/R benchmark images: 

1. Random selection of a real interferometric image 
from the data archive of the MO JAVE prograrrQ; 

2. Selection of a region of this image where no apparent 
sources are present; 

3. Extraction of a noise image of a given size from that 
region; 

4. Construction of a synthetic image composed by N s 
Gaussian sources, with the central one having the 

2 As has been mentioned in Section 2.2, the performance function S* ri 

3 http: //www. physics .purdue . edu/MOJAVE/index . html 



same characteristics as the CLEAN beam parame- 
ters found in the header of the noise image; 

5. Addition of the noise and synthetic images to form 
our new benchmark map. 

Note that this procedure simulates accurately the ther- 
mal noise characteristics of the image, as well as some of 
the CLEAN algorithm effects. The effects of wide gaps 
in the interferometric (u,v) coverage, (e.g., aliasing, si- 
nusoidal ripples, etc.) are not necessarily present in our 
benchmark images. However, a rigorous analysis of real 
CLEAN images can reveal the presence of such (usually 
weak) artifacts (e.g jEkersfl999tlLister et alJl200lD . allow- 
ing the adoption of some additional strategy to avoid mis- 
interpretation of the data image. As regards rippling fea- 
tures, changes in the CLEAN algorithm can substantially 
reduce or even remove the m from the final CLEAN i mage 
(lCornwelllll98l . Similarly. Ipfleiderer k Koebl (fl998l) pro- 
posed a different deconvolution method to translate raw 
complex visibility data into an image that eliminates alias- 
ing effects. 

For the synthetic sources, we decided to use the same 
jet knots Jl and J2 as used in the previous section. The 
noise pattern, and the resulting new benchmark images 
are shown in Figures [5] and HJ Although this procedure 
might not produce a synthetic image with precise repre- 
sentation of the expected noise pattern, we believe that 
it is the best way to generate a more realistic image to 
evaluate the influence of noise in the CE model-fitting op- 
timization process. 

In analogy with the previous section, our CE model- 
fitting algorithm examined the whole image in the x— and 
y— coordinates, searching for Gaussian sources with struc- 
tural parameters in the same range as used in the noiseless 
cases. To guarantee an unbiased solution we ran the CE 
algorithm for different numbers of sources, varying from 
two up to seven. As in the noiseless case, we adopted 
Nq — 30 in the optimization of the Jl image. For J2, 
we had to increase N from 30 to 50 to guarantee a good 
performance of our algorithm. 

Following Section 2.3, we obtained the optimal values for 
the model parameters and their respective uncertainties, 
which are displayed in Tables [5] and [5] As in the noiseless 
cases, comparison between Tables [5] and [5] reveals that our 
algorithm fits tentatively the most intense sources present 
in the image when the number of sources assumed in the 
CE model fitting is smaller than the correct value. If the 
assumed number of the sources is larger than the right one, 
the algorithm puts the extra dim Gaussian sources close 
to the location of the most intense sources. 

As an example, we show in Figure [7] the evolution of the 
values of the CE-optimized model parameters of the noisy 
jet J2 as a function of the iteration number. The param- 
eter values reach a steady-state-like behavior at iterations 
below fc max after converging to the expected values. Simi- 
lar findings are also observed for the other benchmark tests 
discussed in this work. 

The CE model-fitting results are shown in Figures [8] 
and[S] We note that the choice N s = 3 — 6 for jet Jl and 

d has units of Jy 6 beam -6 . 
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N s = 4, 5, and 7 for jet J2 provides more homogeneous 
and less peaked residual maps. 

Considering the uncertainties, there is a plateau after 
A s = 3 in the S^ d graph in the top panel of Figure [TO] 
This means that only a lower limit for N s can be derived 
from such plots. The bottom panel of Figure [TU] points out 
that residuals are minimized for 3 < A s < 6. To break the 
degeneracy in N s , we analyze the values listed in Table [5] 
As in the case of the noiseless images, the optimized pa- 
rameters of the three sources found for N s = 3 also appear 
when we assumed 4 < N s < 6. The extra components 
are dimmer and sometimes superpose completely with the 
brightest component (core). Thus, N s = 3 is fully favored 
in relation to the other possibilities. 

In the case of jet J2, there is also a plateau-like feature 

1 /6 i 

for 7V S > 5 in S ' ' d graph in the top panel of Figure [TT] 
Solutions with N s = 5 and 7 are slightly favored over the 
others in the case of the mean and maximum residuals 
shown in the bottom panel of Figure ITT] The combination 
of both results suggests the existence of five or seven com- 
ponents in jet 32. Comparison between the values listed in 
Table [6] reveals that components 1-5 for N s = 5 were also 
found in the CE optimization using N s — 7 (components 
1, 2, 3, 6, and 7, respectively). As in the case of the noisy 
jet Jl, the extra two components are too dim, compatible 
with null intensity considering the uncertainties. Thus, a 
five-components solution is favored, as expected for this 
benchmark test. 

3.3. Checking the Auto- consistency of the Error 
Estimation of the CE Model-fitting Parameters 

The equations to calculate the model-fitting parameters, 
as well as their respective uncertainties, were introduced in 
Section 2.3. To show that Equations (10) and (11) provide 
a reasonable estimation of the model parameters and their 
errors, we used the three optimization runs performed for 
the noisy jet Jl for A s = 3. The reason for not choosing 
a noiseless image is to keep this analysis as close as possi- 
ble to real situations, which we know to have a structured 
noise pattern embedded into the source signal. 

From the 18,000 best-value data generated in this test 
(three optimization runs, each one composed of 6000 iter- 
ations), we determined the distribution of the logarithm 
of Sp I0 d in the (x, y)-, (e, &)-, and (a, Io)-planes for each of 
the three Gaussian sources, as is displayed in Figure [TSffl . 
Although the S vm & distribution is very irregular, showing 
several local minima in all planes, the CE algorithm was 
able to find the global minimum of S'prod, returning the 
correct parameter values of the ellipses if the error bars 
were considered. For example, the differences between the 
x— and y— values listed in Table [1] and those calculated 
from our technique are much smaller than one pixel in 
both directions. The larger discrepancies were found in 
the (e, 9) — planes but, even in these cases, such differences 
remain within the 3cr errors. 

Our results show also that it is possible to calculate sta- 
tistically the values of the model-fitting parameters and 
their associated uncertainties by applying Equations (10) 
and (11) to the tentative solutions generated in the opti- 
mization process. However, as the real nature of the noise 

4 This representation of <S P rod was chosen due to its complex functional 
test. 



present in interferometric images is not sufficiently known, 
such error estimates should be considered as a lower limit 
for the true uncertainties of the derived Gaussian param- 
eters. Error estimations based on the properties of the 
restoring beam or differences in the fits to t wo or more 
observations acquired ve ry close in time fe.g.. iPiner et al.l 
120071: iLister et alJl2009bh may be used as complements or 
alternatives to Equations (10) and (11) in general situa- 
tions. 

3.4. The Influence of S/R on the CE Model-fitting 
Optimization Process 

The images in the noiseless validation tests had a formal 
infinite S/N, while for the Jl and J2 noise images S/R was 
^100. To check the performance of our CE model-fitting 
algorithm in terms of S/R, we created synthetic images 
composed of the same noise background shown in the top 
panel of Figure [5] and an elliptical Gaussian source cen- 
tered at x = 17.2 pixels and y = 48.1 pixels. We placed 
this source at this particular position in order to super- 
pose it on a region in which there is a strong gradient 
in the noise intensity in the southeast to northwest direc- 
tion, which introduces extra difficulty to the model-fitting 
procedure. The additional parameters were chosen to be 
a = 4.8 pixels, e = 0.85, and 9 = 67?3. The peak intensity 
of this source was varied from 1 to 50 mJy beam -1 , which 
corresponds to S/R ranging from 2 to 100. An example of 
one of these control images is displayed in Figure [TBI 

The efficiency of our algorithm in recovering the ellipti- 
cal Gaussian source parameters as a function of the S/R 
can be seen in Figure [TJ] The values of the Gaussian-peak 
coordinates x and y were recovered with a relative error 
below the 10% level for S/R> 10. The same error level 
was achieved for the parameters e and 9 for S/R> 20, Jo 
for S/R> 30, and a for S/R> 50. For SNR = 100 all pa- 
rameters were found with a relative error below the 5% 
level. 

In terms of the absolute values and their respective un- 
certainties, our CE model fitting technique was able to 
recover the parameter 1$ within the 3cr level for S/R> 20, 
while for the other parameters the same applied for S/R> 
10. 

3.5. Comparing the CE Benchmark Results with Those 
Obtained from the AIPS Task IMFIT 

The aim of this subsection is to compare our findings 
with those obtained from the traditional and widely used 
two-dimensional Gaussian fitting task IMFIT included 
in the Astronomical Image Processing System Package 
(AIPS). 

A comparison between CE and IMFIT model-fitting re- 
sults in the case of the S/R limit tests is shown in Figure 
1141 Taking into account the uncertainties, both techniques 
have very similar behavior in terms of recovering the source 
parameters for this particular test. 

We applied the task IMFIT to the image of the noisy jet 
Jl displayed in Figure [5] A comparison between the val- 
ues of the parameters of the three Gaussian sources of Jl 
from IMFIT and from our CE technique is shown in Table 
[7J Again, both algorithms found the Gaussian parameters 

form, which depends on 18 parameters for this particular benchmark 
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with similar high efficiency. 

Unfortunately, IMFIT cannot fit simultaneously more 
than four Gaussian sources. To make possible a compari- 
son between those algorithms, we eliminate source number 
3 of the noisy jet J2 transforming it into a four-component 
image labeled as JS0. This makes J3 map as a set of two 
overlapping components with S/R ~ 100. We show in Ta- 
ble [8] the model- fitting values recovered from the CE and 
IMFIT fittings as a function of the correct ones. Unlike the 
earlier comparisons presented in this subsection, the CE 
technique had a superior performance in determining the 
structural parameters of the Gaussian sources, especially 
in the case of the dimmest component for which IMFIT 
failed to converge to the correct parameters. In terms of 
relative errors, our CE algorithm was able to recover 7 out 
of 24 parameters within a zero-percent level, and only one 
parameter with an error higher than 5 % level (more pre- 
cisely, the 6.8% level for the parameter 9 of the dimmest 
source). On the other hand, the relative errors obtained 
from IMFIT are systematically higher in 19 out of 24 pa- 
rameters. In the case of the dimmest source, the relative 
errors ranged from 16% to 400%, indicating that the IM- 
FIT did not converge for this component. 

In summary, comparison between the CE optimization 
algorithm and IMFIT suggests that both tools have sim- 
ilar performance in dealing with images that are not too 
complex (few and relatively well-separated sources). How- 
ever, our technique has a better performance when the 
intcrfcromctric map is more complex, as in the case of test 
J3. Note also that our method does not have a limit for 
the number of components to be fitted simultaneously in 
the image, as in the case of IMFIT. 

4. APPLICATION OF THE CE TECHNIQUE TO A REAL 
INTERFEROMETRIC IMAGE 

We have presented in previous sections some valida- 
tion tests necessary to demonstrate the capability of our 
model-fitting technique of dealing with interferometric im- 
ages. The objective of this section is to check the be- 
havior of our algorithm when applied to real interfer- 
ometric images. To do this, we selected a naturally 
weighted /—i mage taken from th e MOJAVE/2cm Survey 
Data Archive ()Lister et al.ll2009al) that corresponds to the 
15 GHz radio map of the BL Lac OJ287 obtained in 1996 
May 27. 

4.1. CE Model-fitting and the Number of Jet Components 

The original fits image of the OJ287 is formed by an 
array of 512 x 512 pixels but only a relatively small frac- 
tion has a jet signal significantly higher than the noise. 
Because of this, we decided to crop the original fits image 
to maintain only the fraction with a useful signal, which 
meant a 66 x 51 pixel image. It is important to empha- 
size that this reduction helps the algorithm to find the 
optimal solutions more rapidly since the parameter space 
is substantially narrowed in this case. 

Several works in the literature have assumed a circular 



Gaussian shape for jet features (e.g., Lobanov et al.|[200lt 
Uorstad et alfl2005t lAgudo et all 120071) . In this work, we 
have assumed a two-dimensional elliptical Gaussian func- 

5 The IMFIT failed to fit four components in the original image of the 
issue. Using three components, IMFIT returns parameters very similar 



tion as being representative of the brightness distribution 
of the jet knots, which will be kept for the radio map of 
the BL Lac O J 287. The justification for this is that our 
CE model fitting will be executed in the image plane, in 
which the jet components are convolved with the elliptical 
synthesized CLEAN beam of the interferometric experi- 
ment. 

We applied our CE algorithm to the cropped image, 
varying the number of elliptical Gaussian sources from two 
to seven. For each adopted number of sources we run the 
algorithm three times. The model fitting processes ex- 
amined the whole image in the x— and y— coordinates, 
searching for Gaussian sources with structural parameters 
in the following ranges: 4.0 < a(pixel) < 30, < e < 0.9, 
—90° < 6 < 90° (positive values for north to east direc- 
tion), and 9.51 x 10~ 4 < I (Jy beam- 1 )^ 1.2. The lower 
value of Iq corresponded to twice the nominal root mean 
square of the O J 287 image. 

The optimal values for the model parameters and their 
respective uncertainties were estimated following Section 
2.3 and are displayed in Table [51 The CE optimization 
results can be seen in Figure 1151 in which the Gaussian 
components are shown superposed on the observed image, 
as well as the residual maps. A visual inspection of the 
residual maps reveals less-structured residuals when three 
to five sources are employed in the optimization. The top 
panel of Figure [TC] shows that the performance function 
introduced in Section 2.2 is better minimized by adopting 
N s > 3. The behavior of the mean and maximum values of 
the residual image shown in the bottom panel of Figure [TBI 
seems to be minimized for 3 < iV s < 5. Therefore, the use 
of Figures [T5] and [TC] only provides an optimal range for 
the number of components present in the image of O J 287. 
However, looking carefully at the entries of Table [51 we 
realize that components 1, 2, and 3 for N s = 3 were also 
found in the CE optimizations using N s = 4 and 5. The 
extra components in both cases are too dim (~ 1 mJy 
beam -1 ), compatible with the zero-intensity level if un- 
certainties are considered. Besides, similar plots of those 
shown in Figure [7] reveal that those extra components in- 
troduced large oscillations in the values of the Gaussian 
parameters. We believe that such findings point firmly to 
the existence of three components in the jet of O J 287 in 
1996 May 27. 

4.2. Comparing the CE Results with Those in the 
Literature 

As discussed in the previous section, our results seem 
to indicate the presence of three components in the jet 
of O J 287 (core pl us two components). Interestingly, 
iLister et all (|2009bl) also assumed the existence of three 
sources in the 15 GHz image of OJ 287 in 1996 May 27. 
Their modeling results were obtained by adopting circular- 
shaped sources fitted in the interferometric (u, v) plane in- 
stead of in the image plane, as has been done in this work. 
A comparison between our mo del-fitting paramete r values 
for N s = 3 and those found bv lLister et all (|2009bf ) can be 
seen in Table [TOl 

To construct Table [TU1 we transformed our original op- 
timized data listed in Table [HI to the format given in Table 

noisy jet J2, displaying an error message related to some convergence 
to those listed in the three-components entry of the Table [6] 
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1 of iLister et al.l (|2009bl ). The first parameter is the core- 
component distance r that can be calculated from 

r = \/Jx- Xcore) 2 + {V ~ 2/core) 2 , (12) 

where x corc and y CO rc are, respectively, the right ascension 
and declination coordinates of component 1 in Table |H1 

The position angles of the components me asured in re- 
lation t o the core component, 77, taken from ILTster et ail 
(2009b) were subtracted by 360° before being displayed in 
Table US 

As mentioned previously, we fitted e ll iptical Gaussian 
sources in the image plane. ILister et al.l (|2009bh adjusted 
circular Gaussian sources to the (it, v) data, which al- 
lowed them to obtain for some components sizes smaller 
than that of the restoring beam. Remembering that the 
convolution of two Gaussian functions results in a Gaus- 
sian function with a squared full width at half-maximum 
(FWHM) proportional to their FWHM, we determined 
the effective size of the optimized components ofwhm 
through: 

"fwhm = \Ja 2 {\- e 2 ) 1/2 - dpWHM ^FWHM ' ( 13 ) 

where 6p^ HM and &^whm are respectively the FWHM of 
the major and minor sy nthesized beam axes for the image 
epoch (1.16 x 0.60 mas: ILister et al.ll2009bh . 

The flux density F of the jet components can be esti- 
mated from: 

F = 81n2 (^viIlM /o . (14) 

\ "fwhm"fwhm / 

In the last equation, a, 6p^ HM and £>fwhm mu st be 
given in pixels, and Iq in Jy beam -1 in order to have F in 
units of Jy. 

We note that there is a relati ve good agreement be tween 
our results and those found by ILister et al.l (|2009bl) if we 
take into account the uncertainties. It suggests that our 
CE technique may provide results as robust as those model 
fittings performed in the (u, v) plane. Further studies us- 
ing a larger set of real images are required to confirm un- 
doubtedly this claim. 

The optimized Gaussian sources shown in Figure [15] 
represent the convolution between the source itself and 
the synthesized beam of the observation. If the source 
is punctual, such a convolution must return the shape of 
the synthesized beam. For the image of OJ 287 consid- 
ered in this work, the beam is elliptical, with its major 
and the minor semi-axes being respectively 4.91 and 2.56 
pixels (or 0.58 and 0.30 mas) in size, and wit h a position 
angle P. A. = 22?28 on the plane of the sky (jLister et al.l 
2009b). This last quantity is related to our definition of 
9 through P.Aa. = 90° + 8. A comparison of these val- 
ues with those listed in Table for N s = 3 shows that 
component 1 is unresolved by the observation, presenting 
major and minor semi-axes equals respectively to 4.92 and 
2.59 pixels, as well as P. A. = 22?70, in full agreement with 
what is expected from a punctual source, which is indeed 
the compact core of OJ 287. This result further reinforces 
the great potential of our technique in modeling interfero- 
metric radio images. 



5. CONCLUSION 

We have developed a new method to obtain model fit- 
tings to interferometric radio images of astrophysical jets 
using a global optimization algorithm known as cross- 
entropy. To validate our model-fitting optimized proce- 
dure, we built benchmark tests that employed synthetic 
images created to simulate as realistically as possible typ- 
ical real interferometric radio maps. We assumed that 
the individual sources of a jet can be represented by two- 
dimensional Gaussian functions defined by six parameters: 
peak intensity coordinates, peak intensity, size and angle 
of the semi-major axis, and eccentricity. 

The first validation tests optimize synthetic images with 
infinite S/R composed of three and five sources (jets Jl and 
J2, respectively). The parameters for each source were se- 
lected in an attempt to reproduce typical image character- 
istics encountered in VLBI AGN jet images. The results of 
our fitting technique for these synthetic images were excel- 
lent: all parameter values were recovered by our technique 
taking into account the formal uncertainties. The relative 
errors of these two benchmark tests are null in most cases, 
and do not exceed 0.25% in the worst cases. 

The second set of benchmark tests was built from jets Jl 
and J2, embedding them in a structured noise backgrounds 
extracted from typical MOJAVE images. Our CE model- 
fitting technique was able to recover the parameters of the 
three sources of the noisy jet Jl in most cases with null 
relative error. In the case of the complex noisy jet J2, 18 
out of 30 parameters had values below of 0.3% in terms of 
relative errors (seven among them having a zero-percent 
value) . 

Plots of the values of the performance function and the 
mean and maximum amplitude of the residual images as a 
function of the number of tentative components only pro- 
vided an optimal range for N s in the majority of the bench- 
mark tests presented in this work. However, we showed 
that a careful analysis of the optimized parameters among 
all used values of N s can remove such degeneracy, and de- 
termine the true value of the number of components. If 
the assumed N s is smaller than the real one, our CE opti- 
mization tries to fit the brightest components since their 
contribution must be preponderant in the minimization of 
the residuals. On the other hand, if the adopted N s is 
larger than the correct one, the extra components tends 
to be too dim (compatible with null intensity if uncertain- 
ties are considered) or they are practically coincident with 
the most intense sources in the image. Thus, the true jet 
components will always be present in the optimized im- 
ages when we assumed a number of sources greater than 
the real one. 

The number of components also has an influence on the 
convergence of the parameters: the inclusion of more com- 
ponents than necessary introduces an oscillatory pattern, 
as well as some strong discontinuities in the plots of the 
values of the Gaussian parameters as a function of the it- 
eration number. This can also be used to constrain the 
value of N s . 

The cross-entropy technique is able to recover the pa- 
rameters of the sources with a similar accuracy to that 
obtained from the traditional AIPS task IMFIT when the 
image is relatively simple (e.g., few components). For more 
complex interferometric maps, our method exhibits a su- 
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perior performance in recovering the parameters of the jet 
components. 

To verify the performance of our algorithm in the case of 
real observational data, we applied it to the image of the 
BL Lac OJ 287 obtained in 1996 May 27 by the MOJAVE 
consortium, as mentioned in Section 4. The behavior of 
the objective function, as well as the mean and the maxi- 
mum values of the residual maps as a function of the num- 
ber of the sources, indicates the presence of three to five 
components. This degeneracy is broken when the data of 
Tableware carefully analyzed. They indicate the presence 
of the three sources found in the optimization with N s = 3 
in the results from N s = 4 and 5, with the extra compo- 
nents too dim to be taken seriously. It is important to em- 
phasize that the structural parameters of the core found 
by our technique are in full agreement with those expected 
from the convolution of a point- like source a n d an el liptical 
restoring beam. Interestingly, iLister et al.l ()2009bf ) mod- 
eled the same image of O J 287 in the (u, v) plane as the 
composition of three circular Gaussian sources. The struc- 
tural source parameters obtained by these authors are in 
good agreement with those found in this work from our 
CE optimization technique. This suggests that our image- 
based model-fitting method might be as efficient as (u, v)- 
based modeling despite the usual deconvolution and bad 
(u, v) coverage artifacts that might be present in CLEAN 
images. The next step in our tests will be to adapt the CE 
method to fittings in the (u, v) plane. 



It is important to point out that we have assumed 
throughout this work that the components can be mod- 
eled by two-dimensional Gaussian functions, which might 
not be suitable in some situations (e.g., in pronounced 
bow-shock regions usually seen in terminal jets) . Indeed, 
there is no prior guarantee that elliptical Gaussian func- 
tions can always represent the correct shape of jet knots 
in general. 

We believe that our results indicate that this new op- 
timization technique can provide a major contribution in 
obtaining model fittings to VLBI images of astrophysical 
jets. In particular, our CE model-fitting technique must 
be used in studies of complex emission regions presenting 
more than three sources, even though it is substantially 
slower than current model-fitting tasks (at 10,000 times 
slower for a single processor, depending on the number of 
sources to be optimized). As in the case of any model 
fitting performed in the image plane, caution is always 
required in analyzing images constructed from a poorly 
sampled (it, v) plane. 

This work was supported by the Brazilian Agencies 
FAPESP (Proc. 2006/57824-1) and CNPq. This research 
has made use of data from the MOJAVE database that is 
maintained by the MOJAVE team (Lister et al., 2009a). 
The authors acknowledge very helpful remarks from an 
anonymous referee. 
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Fig. 1. — Image of the synthetic jet built from three discrete elliptical Gaussian components (see Table[l]for the jet component parameters). 
The gray-scale map as well as the contours is in logarithm scale in units of Jy beam -1 . 
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Fig. 2. — Image of the synthetic jet built from five discrete elliptical Gaussian components (see Table[2]for the jet component parameters). 
The gray-scale map as well as the contours is in logarithm scale in units of Jy beam -1 . 
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Fig. 3. — Quantitative analyses of the model fittings in the case of the noiseless synthetic image of Jl. Top: one-sixth of the performance 
function as a function of the number of sources used in the CE optimization. Taking into account the error bars, there is a plateau-like 
structure for N B > 3. Bottom: Behavior of the mean value of the residuals in terms of the number of sources used in the CE optimization. 
The error bars are the standard deviation of the corresponding residual images. The dotted line refers to the maximum value of the residuals. 
The mean and maximum values reach their minimum at N B = 3. 
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Fig. 4. — Quantitative analyses of the model fittings in the case of the noiseless synthetic image of J2. Top: One-sixth of the performance 
function as a function of the number of the sources used in the CE optimization. Taking into account the error bars, there is a plateau-like 
structure for N a > 5. Bottom: Behavior of the mean value of the residuals in terms of the sources used in the CE optimization. The error 
bars correspond to the standard deviation of the residual images. The dotted line refers to the maximum value of the residuals. 
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Fig. 5. — Left: Noise pattern extracted from a MOJAVE project image used to build the background noise. RMS of this noise image is 
0.515 mjy beam -1 . Right: Image of the noise validation test composed of the summation of the contribution of three elliptical Gaussian 
sources of jet Jl and the background noise shown in the left panel. 
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Fig. 6. — Left: Noise pattern extracted from a MOJAVE project image used to build the background noise. RMS of this noise image 
is 0.463 mjy beam" 1 . Right: Image of the noise validation test composed of the summation of the contribution of five elliptical Gaussian 
sources of jet J2 and the background noise shown in the left panel. 
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Fig. 7. — Evolution of the model parameters of the Gaussian sources of the noisy jet J2 as a function of the iteration number k. The solid 
horizontal lines represent the expected values for that benchmark test (see Table 
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FlG. 8. — Results of the model fitting for the synthetic image of the jet Jl with the addition of a realistic noise pattern obtained after 
varying the number of sources from two to seven (left- to- right direction). Top : Contour lines are of the original image shown in Figure [5] 
the gray-scale image is constructed from the fitted source parameters, and the dark ellipses are the contours of the individual fitted sources 
(respective centers marked with crosses) at the FWHM. Bottom : Respective residual maps in linear scale. 
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FlG. 9. — Results of the model fitting for the synthetic image of jet J2 with addition of a realistic noise pattern obtained after varying the 
number of sources from two to seven (left-to-right direction). Top : Contour lines are of the original image shown in Figure [6] the gray-scale 
image is constructed from the fitted source parameters and the dark ellipses are the contours of the individual fitted sources (respective centers 
marked with crosses) at the FWHM. Bottom : Respective residual maps in linear scale. 
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Fig. 10. — Quantitative analyses of the model fittings in the case of the synthetic noise image of jet Jl. Left: One-sixth of the performance 
function as a function of the number of sources used in the CE optimization. Taking into account the error bars, there is a plateau-like 
structure for N s > 3. Right: Behavior of the mean value of the residual image shown in Figure [8] in terms of the sources used in the CE 
optimization. The error bars correspond to the standard deviation of the residual images. The dotted line refers to the maximum value of 
the residuals. 
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Fig. 11. — Quantitative analyses of the model fittings in the case of the synthetic noise image of jet J2. Left: One-sixth of the performance 
function as a function of the number of sources used in the CE optimization. Right: Behavior of the mean value of the residual image shown 
in Figure [9] in terms of the sources used in the CE optimization. The error bars correspond to the standard deviation of the residual images. 
The dotted line refers to the maximum value of the residuals. 
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Fig. 12. — Mapping of the logarithm distribution of S pr od in terms of the CE elliptical Gaussian parameters for components 3, 2, and 
1 (rightmost, middle, and leftmost panels, respectively), represented by the gray scale (the lower the value of S pro d, the darker the color). 
Square symbols mark the parameter values listed in Table [T] while open circles represent the CE model fitting parameters of the three-source 
entries in Table [5] Error bars correspond to the 3a uncertainties. Note that the CE values are in agreement with the real ones and both are 
located at the deepest minimum of the objective function in each image. 
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Fig. 13. — Image of a synthetic elliptical Gaussian component superposed on the background noise shown in the left panel of Figure [5] with 
S/N equal to 30. The gray-scale map as well as the contours is in linear scale in units of Jy beam -1 . 
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Fig. 14. — Relative errors in the parameters of an elliptical Gaussian source as a function of the S/N. The black circles correspond to the 
results obtained from our CE model fitting technique, while the open squares refer to the model fitting obtained from the AIPS task IMFIT. 
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Fig. 15. — Model fittings in the case of 15 GHz VLBA image of the BL Lac OJ287 using elliptical Gaussian sources. Left : White contour 
lines represent the VLBA map of OJ 287 in 1996 May 27, the gray-scale image is constructed from the fitted source parameters and the dark 
thick ellipses are the contours of the individual fitted sources at the FWHM (center marked with crosses). The number of fitted Gaussian 
sources increases from the top to the bottom panel (from two to seven sources, respectively). The thin white ellipses at the lower leftmost 
corner in each panel are the FWHM synthesized beam of the observation. Right : Residual maps in linear scale. The lower left corner in 
all panels, i.e. the origin of the right ascension and declination offset coordinates, represents pixels 240 and 230 of the original fits image of 
OJ 287 in right ascension and declination directions, respectively. 
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Fig. 16. — Quantitative analyses of the model fittings in the case of 15 GHz VLBA image of the BL Lac OJ 287 in 1996 May 27 using 
elliptical Gaussian sources. Top: One-sixth of the objective function as a function of the number of the sources used in the CE optimization. 
Note the plateau-like structure (considering the error bars) for N B > 3. Bottom: Behavior of the mean value of the residual image shown in 
the lower panels of the Figure fT51 in terms of the sources used in the CE optimization. The error bars correspond to the standard deviation 
of the residual image. The dotted line refers to the maximum value of the residuals. 
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Table 1 

Parameters of the Three Gaussian Sources in the Case of the Jet Jl, Shown in Figure [T] 





X 


y 


a 


e 


e 


Iq 




(pixel) 


(pixel) 


(pixel) 




(dog) 


(Jy/beam) 


1 


11.0 


71.0 


5.5 


0.90 


75.0 


2.0 


2 


17.0 


54.0 


6.5 


0.80 


83.0 


0.3 


3 


13.0 


44.0 


8.0 


0.70 


41.0 


0.2 



Note. — Component 1 represents the core, while components 2 and 3 are the jet knots. 

Table 2 

Parameters of the Five Gaussian Sources in the Case of the Jet 32, Shown in Figure [2j 





X 


V 


a 


e 


9 


Io 




(pixel) 


(pixel) 


(pixel) 




(dog) 


(Jy beam -1 ) 


1 


25.0 


52.0 


3.0 


0.60 


-70.0 


1.00 


2 


27.0 


46.0 


3.2 


0.65 


-50.0 


0.50 


3 


33.0 


38.0 


4.0 


0.70 


-10.0 


0.10 


4 


36.0 


28.0 


5.0 


0.75 


70.0 


0.30 


5 


23.0 


20.0 


8.0 


0.80 


20.0 


0.05 



Note. — Component 1 represents the core, while components 2-5 are the jet knots. 

Table 3 

Model Fitting Parameters of the Gaussian Sources in the Case of the Noiseless Synthetic Image of Jl. 



X 


y 


a 


e 


9 


Io 


(pixel) 


(pixel) 


(pixel) 




(deg) 


(Jy/beam) 



1 


11.00 ± 0.01 


71 


.02 


± 





.02 


5.47 ± 0.02 


0. 


90 


± 





.00 


74.96 ± 0.16 


2. 


.003 


± 





.006 


2 


15.56 ± 0.04 


49 


.55 


± 





.06 


9.21 ± 0.07 


0. 


.85 


± 





.00 


69.39 ± 0.48 


0. 


.368 


± 





.002 


1 


11.00 ± 0.00 


71 


.00 


± 





.00 


5.50 ± 0.00 





90 


± 





.00 


75.00 ± 0.01 


2. 


.000 


± 





.001 


2 


17.00 ± 0.01 


54 


.00 


± 





.02 


6.50 ± 0.01 


0. 


80 


± 





.00 


83.00 ± 0.17 


0. 


.300 


± 





.001 


3 


13.00 ± 0.02 


44 


.00 


± 





.02 


8.00 ± 0.02 


0. 


70 


± 





.00 


40.99 ± 0.39 


0. 


.200 


± 





.001 


1 


11.00 ± 0.00 


71 


.00 


± 





.00 


5.50 ± 0.00 


0. 


90 


± 





.00 


75.00 ± 0.02 


2. 


.000 


± 





.001 


2 


17.00 ± 0.01 


54. 


.00 


± 





.02 


6.50 ± 0.02 


0. 


.80 


± 





.00 


83.00 ± 0.20 


0. 


.299 


± 





.001 


3 


17.54 ± 4.74 


52 


.96 


± 


9 


.80 


5.13 ± 2.73 


0. 


.54 


± 





.13 


40.02 ± 102.94 


0. 


.001 


± 





.001 


1 


13.00 ± 0.03 


44. 


.00 


± 





.04 


8.00 ± 0.03 


0. 


70 


± 





.00 


41.02 ± 0.43 


0. 


.200 


± 





.001 


1 


23.88 ± 0.32 


87 


.79 


± 





.23 


22.97 ± 0.55 





.94 


± 





.01 


89.00 ± 1.54 


0. 


.000 


± 





.001 


2 


26.69 ± 1.57 


86 


.26 


± 


2 


.24 


37.08 ± 2.65 


0. 


.92 


± 





.04 


80.14 ± 10.76 


0. 


.000 


± 





.001 


3 


11.00 ± 0.00 


71 


.00 


± 





.01 


5.50 ± 0.01 





90 


± 





.00 


75.00 ± 0.07 


2. 


.000 


± 





.003 


1 


17.00 ± 0.04 


54. 


.00 


± 





.09 


6.49 ± 0.07 


0. 


.80 


± 





.01 


83.14 ± 0.82 


0. 


.300 


± 





.004 


5 


13.00 ± 0.12 


44. 


.01 


± 





.13 


7.98 ± 0.09 


0. 


70 


± 





.01 


41.10 ± 1.61 


0. 


.200 


± 





.003 


1 


12.93 ± 4.78 


78. 


.87 


± 


9 


.04 


3.41 ± 1.43 


0. 


.63 


± 





.11 


19.44 ± 91.56 


0. 


.001 


± 





.001 


2 


11.00 ± 0.00 


71 


.00 


± 





.00 


5.50 ± 0.00 


0. 


90 


± 





.00 


75.00 ± 0.03 


2. 


.000 


± 





.001 


3 


10.16 ± 3.87 


68 


.13 


± 


6 


.75 


3.44 ± 1.48 


0. 


52 


± 





,19 


30.23 ± 80.18 


0. 


.001 


± 





.001 


1 


9.15 ± 2.44 


63 


.30 


± 


3 


.90 


3.21 ± 1.10 


0. 


61 


± 





.12 


45.78 ± 62.45 


0. 


.001 


± 





.001 


5 


17.00 ± 0.01 


54. 


.00 


± 





.03 


6.50 ± 0.02 


0. 


.80 


± 





.00 


83.00 ± 0.27 


0. 


.300 


± 





.001 


6 


13.00 ± 0.04 


11 


.00 


± 





.04 


8.00 ± 0.03 


0. 


70 


± 





.00 


40.99 ± 0.55 


0. 


.200 


± 





.001 


1 


9.45 ± 0.01 


69 


.86 


± 





.00 


2.81 ± 0.01 





00 


± 





.00 


-6.75 ± 0.29 


0. 


.739 


± 





.005 


2 


3.11 ± 0.05 


65 


.92 


± 





.04 


4.14 ± 0.01 





95 


± 





.00 


51.00 ± 0.11 


0. 


.000 


± 





.001 


3 


12.33 ± 0.02 


64 


.20 


± 





.05 


10.81 ± 0.02 





76 


± 





.00 


57.46 ± 0.22 


0. 


.717 


± 





.003 


1 


1.00 ± 0.06 


56 


.16 


± 





.03 


5.40 ± 0.00 





33 


± 





.00 


-1.58 ± 0.11 


0. 


.343 


± 





.001 


5 


11.62 ± 0.03 


15 


.71 


± 





.02 


29.57 ± 0.09 


0. 


.21 


± 





.00 


90.00 ± 0.33 


0. 


.010 


± 





.000 


6 


1.00 ± 0.03 


38 


.16 


± 





.10 


14.20 ± 0.04 


0. 


79 


± 





.00 


68.25 ± 0.04 


0. 


.000 


± 





.000 


7 


16.70 ± 0.01 


11 


.98 


± 





.11 


2.50 ± 0.02 





16 


± 





.00 


-57.52 ± 0.35 


0. 


.104 


± 





.000 



Note. - 



The number of the sources was varied from two to seven. The uncertainties in each parameter correspond to the 3<r level. 
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Table 4 

Model Fitting Parameters of the Gaussian Sources in the Case of the Noiseless Synthetic Image of 32. 

x y a e Ig 



(pixel) (pixel) (pixel) (deg) (Jy/bcam) 



1 


25.24 


± 


0.05 


52.47 ± 0.07 


2.85 ± 0.11 





.62 


± 





.05 


-72.85 ± 4.00 


0. 


.731 


± 


0. 


.048 


2 


26.18 


± 


0.06 


47.63 ± 0.28 


4.84 ± 0.05 





.86 


± 





.00 


-59.23 ± 1.72 





.579 


± 


0. 


.033 


3 


35.65 


± 


0.06 


28.33 ± 0.08 


6.34 ± 0.08 





.80 


± 





.01 


73.32 ± 1.32 





.279 


± 





.003 


1 


25.19 


± 


0.10 


52.55 ± 0.24 


2.75 ± 0.22 





.57 


± 





.06 


-74.18 ± 5.69 


0. 


.693 


± 





.143 


2 


26.14 


± 


0.39 


47.85 ± 0.96 


4.70 ± 0.46 





.85 


± 





.05 


-61.27 ± 2.85 


0. 


.615 


± 





.080 


3 


35.93 


± 


0.32 


29.44 ± 0.92 


6.32 ± 0.29 





.85 


± 





.03 


-89.87 ± 5.30 


0. 


.234 


± 


0. 


.018 


4 


30.37 


± 


1.45 


23.75 ± 0.85 


10.36 ± 0.73 





95 


± 





.02 


28.27 ± 3.41 





077 


± 


0. 


.024 


1 


25.00 


± 


0.01 


52.00 ± 0.01 


3.00 ± 0.01 





.60 


± 





.00 


-70.03 ± 0.46 


1. 


.000 


± 





.002 


2 


27.00 


± 


0.01 


46.00 ± 0.02 


3.20 ± 0.02 





65 


± 





.01 


-49.98 ± 0.71 





500 


± 





.002 


3 


32.99 


± 


0.09 


38.00 ± 0.07 


4.00 ± 0.07 





.70 


± 





.02 


-10.06 ± 2.67 


0. 


100 


± 





.002 


4 


36.00 


± 


0.02 


28.00 ± 0.04 


5.00 ± 0.03 


0. 


.75 


± 





.00 


70.00 ± 0.60 


0. 


300 


± 


0. 


002 


5 


22.97 


± 


0.24 


19.99 ± 0.15 


7.98 ± 0.21 


0. 


.80 


± 





.02 


20.05 ± 2.33 


0. 


.050 


± 


0. 


.001 


1 


18.68 


± 


2.38 


56.72 ± 20.56 


19.05 ± 6.22 





,94 


± 





.08 


-46.64 ± 37.33 





179 


± 


0. 


072 


2 


24.42 


± 


1.44 


51.32 ± 2.98 


2.52 ± 0.38 





.02 


± 





.35 


52.62 ± 57.48 


0. 


876 


± 


0. 


.210 


3 


10.37 ± 


14.34 


39.41 ± 6.39 


15.86 ± 6.07 





.01 


± 





.40 


5.09 ± 36.28 


0. 


112 


± 


0. 


.105 


4 


38.52 ± 7.56 


27.24 ± 13.86 


2.52 ± 0.27 





.02 


± 





.32 


84.62 ± 86.51 


0. 


333 


± 


0. 


373 


5 


12.94 ± 11.25 


17.34 ± 11.55 


14.63 ± 5.96 





.41 


± 





.16 


29.04 ± 21.81 


0. 


003 


± 


0. 


054 


6 


29.39 ± 0.74 


12.51 ± 14.79 


11.52 ± 4.56 


0. 


.72 


± 





.08 


-48.23 ± 26.84 





155 


± 


0. 


086 


1 


25.00 ± 0.01 


52.00 ± 0.01 


3.00 ± 0.01 





60 


± 





.00 


-70.01 ± 0.56 


1. 


000 


± 


0. 


003 


2 


27.00 ± 0.02 


46.00 ± 0.02 


3.20 ± 0.02 





65 


± 





.01 


-49.98 ± 0.89 


0. 


500 


± 


0. 


003 


3 


32.97 ± 0.12 


38.01 ± 0.08 


3.99 ± 0.09 





.70 


± 





.02 


-10.07 ± 3.30 


0. 


099 


± 


0. 


003 


4 


35.01 ± 6.15 


36.91 ± 7.17 


3.81 ± 2.00 





.59 


± 





.45 


-1.59 ± 93.35 


0. 


001 


± 


0. 


002 


5 


35.86 


± 4.82 


36.75 ± 4.40 


3.57 ± 1.52 





.59 


± 





.42 


-3.48 ± 86.50 


0. 


001 


± 


0. 


002 


6 


36.00 ± 0.03 


27.99 ± 0.04 


5.00 ± 0.04 


0. 


.75 


± 





.01 


69.93 ± 0.73 


0. 


300 


± 


0. 


002 


7 


22.97 ± 0.28 


19.99 ± 0.17 


7.96 ± 0.24 





.80 


± 





.02 


19.99 ± 2.77 


0. 


050 


± 


0. 


001 



Note. - 



The number of the sources was varied from three to seven. The uncertainties in each parameter correspond to the 3<r level. 
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Table 5 

Model Fitting Parameters of the Gaussian Sources in the Case of the Synthetic Image of Jl with Addition of 

Realistic Noise Pattern. 





X 

(pixel) 


y 

(pixel) 


a 

(pixel) 






e 






6 

(deg) 




Io 

(Jy/beam) 


1 


11.00 ± 0.01 


71.04 ± 0.03 


5.44 ± 0.03 





.90 


± 





00 


74.85 ± 0.24 


2 


.005 


± 


0. 


.009 


2 


15.60 ± 0.05 


49.57 ± 0.10 


9.21 ± 0.11 





.85 


± 


0. 


.01 


-86.52 ± 0.82 





.365 


± 





004 


1 


11.00 ± 0.00 


71.00 ± 0.01 


5.50 ± 0.01 





.90 


± 





.00 


75.00 ± 0.08 


2 


.000 


± 





003 


2 


17.00 ± 0.04 


54.02 ± 0.08 


6.48 ± 0.06 





.80 


± 





01 


82.64 ± 0.74 





.300 


± 





.004 


3 


13.02 ± 0.10 


43.97 ± 0.12 


8.02 ± 0.08 





70 


± 





01 


40.89 ± 1.35 





.200 


± 





.003 


1 


11.00 ± 5.85 


71.00 ± 3.17 


5.50 ± 3.22 





.90 


± 


0. 


41 


75.00 ± 44.09 


1 


.382 


± 





.005 


2 


11.00 ± 3.52 


71.00 ± 6.22 


5.50 ± 2.79 





.90 


± 





15 


75.00 ± 30.45 





.110 


± 





.005 


3 


17.01 ± 0.06 


54.06 ± 0.13 


6.48 ± 0.10 





.80 


± 


0. 


01 


82.91 ± 1.16 





.298 


± 





.007 


4 


13.06 ± 0.15 


44.02 ± 0.18 


8.04 ± 0.13 





.70 


± 


0. 


01 


41.63 ± 2.08 





.201 


± 





.005 


1 


11.00 ± 14.89 


71.02 ± 12.04 


5.56 ± 3.08 





.90 


± 





09 


74.94 ± 15.47 


1 


.153 


± 


0. 


.056 


2 


11.04 ± 15.42 


71.01 ± 13.52 


5.33 ± 21.83 





.89 


± 





09 


75.14 ± 96.36 





.017 


± 


0. 


.063 


3 


11.01 ± 11.51 


70.93 ± 12.11 


5.28 ± 3.08 





.88 


± 





.53 


75.29 ± 105.27 





.022 


± 


0. 


.036 


4 


13.10 ± 0.12 


65.60 ± 0.25 


5.74 ± 0.19 





.89 


± 


0. 


02 


76.01 ± 2.27 





.305 


± 





.090 


5 


14.84 ± 0.28 


47.17 ± 0.35 


8.55 ± 0.31 





.81 


± 


0. 


03 


63.56 ± 4.44 





.273 


± 





.010 


1 


19.62 ± 11.54 


74.11 ± 41.60 


3.46 ± 43.82 





.91 


± 


0. 


.54 


73.29 ± 89.52 





.002 


± 





.010 


2 


11.00 ± 20.64 


71.01 ± 26.13 


5.49 ± 34.44 





.90 


± 





17 


74.95 ± 38.78 





.143 


± 





.010 


3 


11.00 ± 2.63 


71.00 ± 3.50 


5.50 ± 1.26 





.90 


± 





.42 


75.00 ± 83.42 


1 


.933 


± 





Oil 


4 


21.01 ± 11.65 


60.90 ± 40.89 


4.81 ± 46.18 





.90 


± 





57 


75.55 ± 91.04 





.001 


± 





.014 


5 


17.00 ± 0.16 


54.01 ± 0.39 


6.47 ± 0.21 





.80 


± 


0. 


02 


-2.09 ± 1.73 





.297 


± 





.017 


6 


13.05 ± 0.42 


44.04 ± 0.60 


8.02 ± 0.29 





70 


± 


0. 


03 


41.73 ± 5.26 





.202 


± 





013 


1 


23.99 ± 4.21 


82.17 ± 5.84 


11.71 ± 2.36 





.90 


± 





34 


77.29 ± 47.42 





.001 


± 





.013 


2 


11.00 ± 1.53 


71.01 ± 2.94 


5.50 ± 2.31 





.90 


± 





13 


74.97 ± 20.25 





.097 


± 


0. 


.010 


3 


11.00 ± 3.67 


71.00 ± 9.60 


5.49 ± 2.05 





.90 


± 





.53 


75.02 ± 80.05 





.961 


± 


0. 


.017 


4 


20.03 ± 5.86 


60.09 ± 15.96 


8.91 ± 3.75 





.59 


± 





.73 


17.44 ± 123.92 





.001 


± 


0. 


.015 


5 


17.03 ± 4.84 


54.14 ± 9.93 


6.55 ± 3.51 





.80 


± 





.49 


82.69 ± 67.82 





.003 


± 





018 


6 


16.94 ± 1.77 


53.87 ± 4.43 


6.32 ± 1.63 





.79 


± 





.29 


-77.18 ± 25.46 





.209 


± 





.019 


7 


13.09 ± 0.42 


44.18 ± 0.48 


7.97 ± 0.35 





.70 


± 


0. 


04 


43.44 ± 6.16 





.208 


± 





013 



Note. - 



The number of the sources was varied from two to seven. The uncertainties in each parameter correspond to the 3(7 level. 
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Table 6 

Model Fitting Parameters of the Gaussian Sources in the Case of the Synthetic Image of J2 with Addition of 

Realistic Noise Pattern. 





X 

(pixel) 


ti 
y 

(pixe 


1) 


a 

(pixel) 












(deg) 




Io 

(Jy/beam) 


1 


25.20 ± 0.05 


52.55 ± 


0.09 


2.80 ± 0.14 


0. 


59 


± 





.07 


-70.28 ± 6.37 





.740 


± 





.084 


2 


26.20 ± 0.12 


47.59 ± 


0.47 


4.62 ± 0.12 





85 


± 





.01 


-59.15 ± 2.57 





.597 


± 


0. 


.047 


3 


35.72 ± 0.04 


28.35 ± 


0.06 


5.93 ± 0.06 


0. 


.77 


± 





.01 


-90.00 ± 0.27 





.277 


± 


0. 


.003 


1 


25.19 ± 0.09 


52.56 ± 


0.20 


2.74 ± 0.18 





56 


± 





.05 


-74.65 ± 5.19 





.685 


± 





.123 


2 


or 1 9 + 33 


47.90 ± 


0.83 


4.72 ± 0.41 


0. 


85 




o 


.05 


_«i ^2 + 2 16 





.620 


± 





.070 


3 


35.93 ± 0.52 


29.46 ± 


1.11 


6.33 ± 0.29 


0. 


85 


± 





.01 


-89.90 ± 2.84 





.234 


± 





.014 


4 


30.35 ± 1.19 


23.72 ± 


0.71 


10.36 ± 0.54 


0. 


95 


± 





.01 


28.21 ± 2.95 





077 


± 


0. 


.018 


1 


24.99 ± 0.01 


52.08 ± 


0.06 


2.97 ± 0.02 


0. 


59 


± 





.01 


-69.26 ± 1.02 





.981 


± 





.016 


2 


26.92 ± 0.06 


46.14 ± 


0.11 


3.32 ± 0.10 





68 


± 





.02 


-52.07 ± 1.24 





514 


± 


0. 


010 


3 


33.09 ± 0.11 


37.92 ± 


0.09 


3.90 ± 0.11 





69 


± 





.03 


-6.34 ± 5.58 





.100 


± 


0. 


002 


4 


36.00 ± 0.03 


28.03 ± 


0.07 


5.01 ± 0.04 


0. 


.75 


± 





.01 


70.07 ± 0.91 





300 


± 


0. 


002 


5 


22.95 ± 0.36 


19.98 ± 


0.20 


8.00 ± 0.34 





.80 


± 





.02 


18.85 ± 2.67 


0. 


.050 


± 


0. 


002 


1 


22.03 ± 4.07 


57.14 ± 


9.16 


2.85 ± 0.44 


0. 


86 


± 





.14 


28.63 ± 69.10 





.039 


± 





.412 


2 


25.53 ± 7.54 


50.15 ± 


10.96 


2.67 ± 1.84 


0. 


06 


± 





.61 


-47.56 ± 22.53 


1 


.120 


± 


0. 


.852 


3 


26.01 ± 7.48 


43.48 ± 15.16 


17.36 ± 11.51 





.04 


± 





.40 


-53.53 ± 30.24 


0. 


.015 


± 


0. 


154 


4 


26.16 ± 1.00 


38.76 ± 11.33 


29.73 ± 22.66 


0. 


.71 


± 





.11 


-88.43 ± 16.76 





.071 


± 


0. 


756 


5 


31.65 ± 4.64 


37.26 ± 


6.21 


20.23 ± 14.37 


0. 


89 


± 





.63 


77.80 ± 130.15 





.014 


± 


0. 


.146 


6 


32.13 ± 6.73 


12.29 ± 


6.60 


2.95 ± 4.79 


0. 


.72 


± 





.10 


-16.67 ± 26.55 





.403 


± 


0. 


296 


1 


45.84 ± 2.79 


55.77 ± 


2.78 


28.72 ± 2.31 


0. 


78 


± 





.25 


-43.40 ± 58.53 


0. 


.000 


± 


0. 


000 


2 


41.63 ± 10.18 


53.99 ± 


5.73 


10.71 ± 23.06 


0. 


.77 


± 





.28 


-22.08 ± 73.25 





.000 


± 


0. 


001 


3 


25.00 ± 0.01 


52.00 ± 


0.01 


3.00 ± 0.01 


0. 


60 


± 





.00 


-70.09 ± 0.63 


0. 


.999 


± 





.003 


4 


27.00 ± 0.02 


46.00 ± 


0.02 


3.20 ± 0.02 


0. 


65 


± 





.01 


-50.32 ± 0.99 





500 


± 


0. 


003 


5 


33.03 ± 0.12 


37.97 ± 


0.09 


3.99 ± 0.10 





70 


± 





.02 


-10.92 ± 3.53 


0. 


.101 


± 


0. 


002 


6 


35.99 ± 0.03 


27.99 ± 


0.05 


4.99 ± 0.04 


0. 


75 


± 





.01 


69.81 ± 0.82 





.301 


± 


0. 


002 


7 


22.85 ± 0.30 


19.93 ± 


0.18 


7.87 ± 0.26 


0. 


79 


± 





.02 


18.48 ± 2.86 





050 


± 


0. 


002 



Note. — The number of the sources was varied from two to seven. The uncertainties in each parameter correspond to the 3(7 level. 



Table 7 

Comparison between the Relative Errors in the Parameters of the Gaussian Sources of the Noisy Jet Jl 
Obtained from our CE Technique and the AIPS Task IMFIT (Values in Parentheses). 





e x 




e a 




£0 






(%) 


(%) 


(%) 


(%) 


(%) 


(%?) 


1 


0.000 ± 0.000 


0.000 ± 0.000 


0.000 ± 0.000 


0.000 ± 0.000 


0.000 ± 0.053 


0.000 ± 0.050 




(0.005 ± 0.013) 


(0.000 ± 0.013) 


(0.012 ± 0.203) 


(0.001 ± 0.057) 


(0.004 ± 0.100) 


(0.000 ± 0.172) 


2 


0.000 ± 0.118 


0.037 ± 0.056 


0.308 ± 0.462 


0.000 ± 0.000 


0.470 ± 0.410 


0.000 ± 0.667 




(0.035 ± 0.268) 


(0.050 ± 0.137) 


(0.303 ± 1.351) 


(0.290 ± 0.928) 


(0.357 ± 1.065) 


(0.120 ± 1.150) 


3 


0.154 ± 0.308 


0.068 ± 0.091 


0.250 ± 0.375 


0.000 ± 0.000 


0.341 ± 1.415 


0.000 ± 0.500 




(0.177 ± 0.942) 


(0.036 ± 0.266) 


(0.209 ± 2.031) 


(0.182 ± 2.520) 


(0.012 ± 4.929) 


(0.165 ± 1.725) 



Note. — The uncertainties in each parameter correspond to the 3<x level. 



VLBI images and cross-entropy optimizer 



23 



Table 8 

Comparison between the Relative Errors in the Parameters of the Gaussian Sources of the Noisy Jet J3 
Obtained from our CE Technique and the AIPS Task IMFIT (Values in Parentheses). 





(%) 


(%) 


Co 

(%) 


(%) 


ee 

(%) 


<$ 


1 


0.000 ± 0.013 


0.058 ± 0.013 


0.333 ± 0.111 


0.000 ± 0.000 


0.543 ± 0.238 


0.700 ± 0.167 




(4.000 ± 0.009) 


(1.941 ± 0.009) 


(0.135 ± 0.187) 


(0.194 ± 0.398) 


(0.354 ± 0.407) 


(0.180 ± 0.159) 


2 


0.148 ± 0.037 


0.130 ± 0.029 


1.250 ± 0.312 


1.538 ± 0.513 


1.420 ± 0.413 


1.200 ± 0.267 




(4.030 ± 0.064) 


(1.735 ± 0.039) 


(2.514 ± 0.727) 


(5.344 ± 1.392) 


(0.604 ± 2.138) 


(49.778 ± 0.318) 


3 


0.000 ± 0.019 


0.071 ± 0.036 


0.000 ± 0.200 


0.000 ± 0.000 


0.100 ± 0.248 


0.000 ± 0.222 




(3.669 ± 0.059) 


(4.843 ± 0.098) 


(8.331 ± 0.685) 


(2.225 ± 0.545) 


(8.801 ± 0.679) 


(1.293 ± 0.530) 


4 


0.043 ± 0.333 


0.050 ± 0.200 


0.500 ± 0.875 


1.250 ± 0.417 


6.800 ± 2.867 


0.000 ± 0.667 




(13.561 ± 0.077) 


(123.795 ± 0.093) 


(59.578 ± 0.304) 


(16.249 ± 0.927) 


(350.870 ± 4.270) 


(398.780 ± 3.180) 



Note. — The uncertainties in each parameter correspond to the la level. 



Table 9 

Model Fitting Parameters of the Elliptical Gaussian Sources in the Case of 15 GHz VLBA Image of the BL 

Lac OJ 287 Obtained in 1996 May 27. 

x y a e 6 Iq 



(pixel) (pixel) (pixel) (deg) (Jy/bcam) 



1 


253.65 ± 0.02 


257.75 ± 0.01 


4.94 


± 


0.01 





85 


± 





.00 


-66.65 ± 0.21 


1 


023 


± 





.021 


2 


257.79 ± 0.26 


257.39 ± 0.03 


5.35 


± 


0.03 





65 


± 





.01 


-43.55 ± 0.93 





.251 


± 





.018 


1 


253.46 ± 0.03 


257.81 ± 0.03 


4.92 


± 


0.01 





85 


± 





.00 


-67.30 ± 0.09 





963 


± 





.032 


2 


257.08 ± 0.39 


257.34 ± 0.04 


5.10 


± 


0.06 





69 


± 





.00 


-56.31 ± 1.91 


0. 


330 


± 





.022 


3 


267.21 ± 1.71 


257.19 ± 0.45 


4.97 


± 


0.18 





67 


± 





.33 


-66.06 ± 19.54 


0. 


.026 


± 





.005 


1 


253.47 ± 0.02 


257.81 ± 0.02 


4.92 


± 


0.01 


0. 


85 


± 





.00 


-67.26 ± 0.08 





.979 


± 





.016 


2 


257.26 ± 0.19 


257.33 ± 0.02 


5.10 


± 


0.03 


0. 


69 


± 





.00 


-56.36 ± 1.09 


0. 


.319 


± 





Oil 


3 


267.64 ± 0.62 


257.22 ± 0.37 


4.94 


± 


0.24 





.74 


± 





.14 


-69.89 ± 8.21 





.024 


± 





.002 


4 


286.83 ± 27.10 


247.32 ± 18.43 


8.24 


± 


3.91 





.34 


± 





.47 


-26.31 ± 101.36 





.001 


± 





.001 


1 


253.48 ± 0.04 


257.79 ± 0.06 


4.92 


± 


0.01 





.85 


± 





.00 


-67.19 ± 0.22 





.993 


± 





.048 


2 


257.41 ± 0.52 


257.36 ± 0.08 


5.11 


± 


0.04 





70 


± 





.03 


-56.06 ± 1.65 





.309 


± 





.030 


3 


267.80 ± 0.49 


257.22 ± 0.65 


4.94 


± 


0.41 





78 


± 





.13 


-69.20 ± 7.94 





.023 


± 





.003 


4 


270.60 ± 13.21 


261.52 ± 8.71 


6.22 


± 


4.06 





60 


± 





.71 


46.32 ± 103.67 


0. 


.001 


± 





.002 


5 


286.76 ± 46.74 


251.15 ± 16.89 


11.51 


± 


16.73 





63 


± 





.78 


4.61 ± 95.79 





.001 


± 





.002 


1 


249.58 ± 35.10 


262.14 ± 5.32 


10.46 


± 


12.18 





40 


± 





.27 


-49.96 ± 171.24 





.209 


± 





.459 


2 


250.44 ± 6.69 


240.54 ± 38.20 


4.49 


± 


0.96 





70 


± 





.34 


-68.66 ± 3.07 


0. 


.343 


± 


1 


.392 


3 


256.97 ± 23.46 


240.91 ± 35.56 


4.28 


± 


1.16 





66 


± 





.20 


-32.52 ± 79.82 


0. 


.121 


± 





.218 


4 


258.72 ± 18.19 


257.04 ± 5.17 


4.52 


± 


2.42 





14 


± 





.29 


-7.00 ± 63.81 





.469 


± 


1 


.033 


5 


265.66 ± 37.73 


261.34 ± 27.97 


6.89 


± 


2.03 


0, 


74 


± 





.90 


43.50 ± 140.40 





.001 


± 





.001 


6 


272.93 ± 34.85 


255.07 ± 5.02 


13.00 


± 


17.47 





39 


± 





.66 


-12.24 ± 97.00 


0. 


.114 


± 





.461 


1 


246.40 ± 1.61 


262.42 ± 1.05 


6.48 


± 


0.35 





.42 


± 





.10 


12.27 ± 17.84 


0. 


.428 


± 





.146 


2 


252.84 ± 2.89 


231.15 ± 6.23 


11.85 ± 1.37 


0. 


.80 


± 





.07 


-23.07 ± 14.37 





.151 


± 





.033 


3 


256.26 ± 2.03 


257.95 ± 0.68 


4.01 ± 0.58 


0. 


21 


± 





.10 


17.13 ± 9.89 





.685 


± 





.153 


4 


264.70 ± 1.48 


263.77 ± 1.50 


13.26 ± 1.63 


0. 


00 


± 





.14 


75.57 ± 14.68 





.001 


± 





.001 


5 


268.39 ± 2.66 


271.38 ± 3.60 


4.02 ± 1.17 





27 


± 





.09 


-0.70 ± 9.55 





.001 


± 





.000 


6 


300.81 ± 9.59 


261.40 ± 0.90 


19.62 ± 3.29 





83 


± 





02 


-17.38 ± 10.27 





.156 


± 





.023 


7 


303.79 ± 8.29 


234.27 ± 5.07 


29.86 


± 5.68 





10 


± 





.12 


69.01 ± 29.72 


0. 


.001 


± 





.005 



Note. — The number of sources was varied from two to seven. Component 1 refers to the compact core of OJ 287. The uncertainties in each 
parameter correspond to the 3<x level. 



24 



Caproni et al. 



Table 10 

Compariso n betwe en the Values of the Parameters Found in This Work for 7V S = 3 and Those Published in 

iLlSTER ET AL.I (|2009bT ) (iN PARENTHESES) FOR THE 15 GHz INTERFEROMETRIC MAP OF OJ 287 OBTAINED IN 1996 MAY 27. 





F 


r 


V 


"FWHM 




(Jy) 


(mas) 


(deg) 


(%) 


1 


0.96 ± 0.03 






0.06 ± 0.00 




(1.11 ± 0.06) 






(0.08 ± 0.01) 


2 


0.49 ± 0.03 


0.37 ± 0.05 


-97.40 ± 1.45 


0.59 ± 0.01 




(0.31 ± 0.02) 


(0.47 ± 0.05) 


(-95.50 ± 9.55) 


(0.33 ± 0.03) 


3 


0.04 ± 0.02 


1.38 ± 0.24 


-92.60 ± 2.66 


0.58 ± 0.19 




(0.07 ± 0.00) 


(1.10 ± 0.05) 


(-92.20 ± 9.22) 


(0.53 ± 0.05) 



Note. — The uncertainties concerning CE results correspond to lcr-value, which was obtained from standard propagation of the errors. 



